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Abstract: We clarify the relationship between the current formalism developed by 

Gyulassy, Kaufmann and Wilson and the Wigner function formulation suggested by 

O. 

cn . Pratt for the 2-particle correlator in Hanbury-Brown Twiss interferometry. When 



0^ 



^ 



applied to a hydrodynamical description of the source with a sharp freeze-out hyper- 
surface, our results remove a slight error in the prescription given by Makhlin and 



> 

in 

o 

r-^ ' Sinyukov which has led to confusion in the literature. 

O' 

Qh! It is widely accepted that if the nuclear matter created in ultra-relativistic heavy- 

Q \ ion collisions attains a high enough energy density, it will undergo a phase transition 



into a quark-gluon plasma. For this reason, it is of great interest to determine the 



^ . energy densities actually attained in these collisions. The total interaction energy of 

H ; 

. 5t , a given reaction can be directly measured by particle calorimeters and spectrometers. 

Although there is no analogous direct measurement for the size of the reaction region, 
Hanbury-Brown Twiss interferometry [|1| provides an indirect measurement in terms 
of the correlations between produced particles. 

Ten years ago, Pratt |Q used the covariant current formulation of Gyulassy, Kauf- 
mann and Wilson |^ to show that the correlations between two particles could be ex- 
pressed in terms of one-particle pseudo- Wigner functions. Although Pratt's derivation 
was non-relativistic, it provided a valuable link between the experimental data and 
many semi-classical event generators whose output came in the form of one-particle 
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distributions. Since that time, different methods have been used to relativistically 
generahze Pratt's result |^, ^ 0, but to our knowledge, the simplest generalization 
(using the covariant current formalism covariantly) has never been published. The 
aim of this letter is twofold: (1) to fill the above void in the literature, and (2) to 
show that by applying the final result to hydrodynamical models with 3-dimensional 
freezeout hypersurfaces, a dispute in the literature about the correct form of the 
2-particle correlator in these models can be resolved. 

The covariant single- and two-particle distributions for bosons are defined by 

Piip) = E^ = E{a^p)a{p)), (1) 

dN 
A(Pa,P6) = EaEi,— — — — = E„Efe(a+(pa)a+(pfo)a(p6)a(pa)) , (2) 

where a'*'(p) (a(p)) creates (destroys) a particle with momentum p. The two particle 
correlation function is then given by 

^^P-P^^-(iV(iV-l))P,(p.)P,(p,)- ^^^ 

Using the classical covariant current formalism of 0, we will show that for a 
general class of chaotic current ensembles the two particle distribution for bosons 
obeys a Wick theorem: 

i^2(Pa,P.) = ^^^^^~ ^^^ (Pi(pJPi(p,) + |5(p„p,)p) , (4) 

where we define the following covariant quantity 



SijPa, Pb) = \lEaEb (a+(pa)a(p6)) . (5) 

We will then show that S is equal to the Fourier transform of a kind of Wigner 
function: 

5(p„ p,) = 5(q, K) = I d'^x e'^^-- Six, K) , (6) 

where the off-shell 4- vector K = j{pa + Pb) is the average of two on-shell (p° = Ei) 
4-momenta, and q = pa — Pb is the off-shell difference of the same two momenta so 



that their scalar product vanishes, K^^q^ = 0. For the special case of p^ = pb, K = pa 
becomes on-shell and 

S{pa,Pa) = S{0,K) = P^{p,) (7) 

It should be noted that eqn. @ involves a 4- dimensional Wigner transform, in contrast 
to the 3-dimensional expression suggested by Pratt [0] which neglects retardation 
effects. 

In [|] it was shown that a classical source current J(x) generates free outgoing 
pions in a state which satisfies 

a(p)|J) = zJ(p)|J), (8) 

where 

J(p) = / , exp[t{Ept - p-x)] J{x) (9) 



is the on-shell Fourier transform of the source J{x), and {J\J) = J d^p\J{p)\'^ = 1- 
For classical currents, the ensemble expectation values in eqns. (P, (0), and (H) can 
then be defined in terms of a density operator p involving the state \J) such that 

{d) = tT{pd). 

Generalizing the result of in order to allow for arbitrary x — p correlations, we 
consider an ensemble of chaotic source currents at positions Xi with momenta pi, 

N 

J(x) = J2 e'^' e-*P'("-"^) Mx - Xi) , (10) 

where 0j is a random phase. The momenta Pi of the sources can, but need not be on 
the boson mass-shell; for example, the source could be a decaying A-resonance with 
3-momentum pj. The on-shell Fourier transform of (|TUp is 



Jip)=J2e^'^^e^^-^^Mp-p^), (11) 

where 

d'^x 



Up -Pr)= I e'^^-^^>^ Mx) (12) 

(27r)32Ep 



is the Fourier transform of Jo{x), and p is on-shell while pi may be off-shelL 
We then choose a density operator such that 

tripO) =Y.PnII d^x,d%p{x,,p,) / -^{J\0\J) (13) 

Ar=0 1=1'' •'^ ^^ 

where p{xi,pi) is the covariant probabihty density of the source points {xi,pi) in phase 
space, and P/v is the probabihty distribution for the number of sources in the reaction. 
These probabihties are normahzed as follows: 

/oo 
d^xd^pp{x,p) = l, Y.Pn = 1. (14) 

N=0 

Using (P) and the above definitions, it is straightforward to show that 
Pi(p)=Ep(|J(p)|2) = {N)EpJd^x,d'p,p{x,,Pi)\Mp-PiW 



{N)EpJd^prp{pi)\Mp-Pi)\'- (15) 



The single particle spectrum is thus obtained by folding the momentum spectrum 
|Jo(p)p of the individual source currents Jo with the 4- momentum distribution of the 
sources, p{p) = J d'^x p{x,p). 

Similarly, if one neglects cases in which two particles are emitted from exactly the 
same point H, one finds: 



A(Pa,Pfc)= ^^^^^^^^^^ EaE, 



{\J{Pa)\'){\J{p,)\') + {r{Pa)J{p,)){r{p,)J{Pa)) 

(16) 
which proves eqn.(§) by way of (H). 

Using eqn.(P), we find the following relationship: 

J*{Pa) J{Pb) = / TTT^TTTTWW exp(-ipa-a;i + ip^-xa) J*(xi) J(x2) 
J [2iiY2^/EaEb 

= I {2nrZkE, ""P^-^^-^ - ^K.y)r{x + \y)J{x - \y) , (17) 



where x = |(xi + X2) and y = Xi — X2- The above relation proves eqn.(|^) as long as 
the following expression for the Wigner function is used: 

S{x, K)=j ^^ e-^^-y {J\x + \y)J{x - \y)) . (18) 

The average on the r. h. s. is defined in the sense of eqn.(|T3|) and can be evaluated 



with the help of the definition (10) to yield 

S{x,K) = {N) fd^zd\p{x-z,q)So{z,K-q), (19) 

where 

So{x,p) = I ^^ e-^^-yj*{x + \y)Ux - \y) (20) 

is the Wigner function associated with an individual source Jq. Thus the one- and 
two-particle spectra can be constructed from a Wigner function which is obtained 
by folding the Wigner function for an individual boson source Jo with the Wigner 



distribution p of the sources. Eqn. (|T9]) is useful for the calculation of quantum statis- 
tical correlations from classical Monte Carlo event generators for heavy-ion collisions: 
{N)p{x,p) can be considered as the distribution of the classical phase-space coordi- 
nates of the boson emitters (decaying resonances or 2-body collision systems), and 
So{x,p) as the Wigner function of the free bosons emitted at these points. Replacing 
the former by a sum of (5-functions describing the space-time locations of the last 
interactions and the boson momenta just afterwards, and the latter by a product 
of two Gaussians with momentum spread Ap and coordinate spread Ax such that 
AxAp > h/2, we recover the expressions derived in 0. 
Using eqns. (^ to (|^), our final result is then 

C(p„p,) = l + i?(q,K) (21) 

where the "correlator" R is given by 

fi(q.K)= J^fa'^)!' . (22) 

S(0,p.)S(0,p») 

Equation (P) is the starting point for a practical evaluation of the above correla- 
tor. It should be noted that due to the on-shell condition of (P), it is impossible to 



reconstruct S{x, K) from the correlator in a model independent way. Thus any ana- 
lysis of data on -R(q, K) necessarily involves suitable model assumptions for S{x, K), 
in particular for the x — K correlations in the source distribution. In most practical 
applications one takes for S{x, K) a classical (on-shell) phase-space distribution. In 
hydrodynamical models, for example, this phase-space distribution is taken as a local 
equilibrium Bose-Einstein distribution localized on a 3-dimensional freeze-out hyper- 
surface S(x) which separates the thermalized interior of an expanding fireball from 
the free-streaming particles on its exterior p|: 

q (r K\ = ^"" + ^ / K^dh,{x')5^''\x-x') 

"^ ' ^ {2'nyi ke^v{l3{x')[K-u{x')-fio.{x')]}-l' ^ ' 

Here Sa and /i^ denote the spin and chemical potential of the emitted particle species 
a, while Ujy(x), /3(x), and (fa^{x) denote the local hydrodynamic flow velocity, inverse 
temperature, and normal-pointing freeze-out hypersurface element. Inserting this 
equation into (0), one obtains the Cooper- Frye formula 

5(0,p) =Pi(p) = f p^d\{x)f{x,p) (24) 

J 2-1 

where we define the distribution function (for clarity we drop the index a for the 
particle species) 

f{x,p) = ^2^^3 exp{/5(a;)[p-M(a;) - fi{x)]} - 1 " ^^^^ 

For the numerator of the correlator, 

\S{q, K)|2 = ^ K^^d^ix) K^d'aM fi^, K) f{y, K) exp[^g-(x - y)] , (26) 



we find an expression which is very similar to the one given in |TD[. There, however, 
each of the two distribution functions under the integral featured on-shell arguments 
Pa and phi respectively, instead of the common (off-shell) average argument K as in 



(PGl). This error in |I0[ can be traced back to an inaccurate transition from finite 
discrete volumes along the freeze-out surface S to the continuum limit [|ll|]. Taking 
over this inaccuracy produces (in particular for very rapidly expanding sources) un- 



physical [T^ oscillations of the correlator around zero at large values of q [T^, |T^ 



which are inconsistent with the manifestly positive definite nature of the correlator 
The symmetric form (|26| ) (in contrast to the asymmetric one given in |T0[) allows 



one to replace the exponential by the cosine and to split the expression into two real 
3-dimensional integrals: 

|^(q, K)|2 = (5i(q, K))^ + (^^(q, K))^ , (27) 

where 



~S,M. K) = /^ K^d\{x) fix, K) I "^^|^;;j;| I . (28) 



cos(g-x) 
sin 

This facilitates the numerical evaluation of the correlator. 
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